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SECTION  I 


INTRODUCTION 

This  report  presents  a  first  order  mathematical  model  of  the  interior  ballistics  of  two-stage 
light  gas  guns.  The  model  and  experimental  data  ore  presented  for  a  helium  system,  hut  con¬ 
ceptually,  any  driver  gas  can  be  used  if  the  thermodynamic  properties  are  known. 

Section  II  describes  the  models  of  the  combustion  chamber,  the  pump  tube,  and  the 
launch  tube,  as  well  as  the  mechanism  of  coupling  between  the  chambers.  The  results  of  the 
analysis  (Equations  11,  20,  and  24)  are  ir.  an  incremental  form  specifically  tailored  for  digital 
machine  computation. 

The  computer  program  generated  from  the  analysis  is  listed  in  Section  III,  along  with 
samples  of  input  and  output  data,  as  well  as  a  test  case. 

Section  IV  illustrates  the  comparison  of  experimental  to  analytic  data  fit  and  discusses 
the  particular  experimental  system  utilized  to  verify  the  generated  mathematical  analysis. 

The  analytic  approach  used  is  essentially  heuristic.  During  this  effort,  engineering  simplicity 
has  been  selected  over  mathematical  elegance.  The  vague  or  undefined  lumped  parameter 
constant  approach  has  been  avoided  as  much  as  possible. 
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SECTION  |i 


ANALYSIS 


1.  SYSTEM  DESCRIPTION 

The  light  gas  system  used  to  generate  and  calibrate  this  mathematical  model  is  shown 
schematically  in  Figure  1  This  system  is  a  typical  two  stage  light  gas  gun  which  uses  helium 
as  t(v»  working  fluid.  The  mechanism  of  operation  is  to  burn  gun  propellant  in  the  combustion 
chamber  until  a  pressure  is  generated  (around  900  psi)  which,  for  this  device,  will  shear  the 
restraining  ring  on  the  piston  and  allow  it  to  travel  into  and  compress  the  helium  gas  on  the 
pump  tube  staqe.  The  gas  being  compressed  eventually  reaches  a  much  higher  pressure  than 
the  driver  gas,  due  to  the  inertia  of  the  relatively  heavy  piston  traveling  at  a  velocity  of 
approximately  2,000  ft/sec.  This  pump  tube  gas  is  compressed  until  it  reaches  a  pressure 
adequate  to  shear  the  restraining  mechanism  on  the  propelled  payload.  The  projectile  is  then 
accelerated  at  high  velocity  down  the  evacuated  launch  tube,  utilizing  the  very  low  pressure 
gradient  decrement  associated  with  the  low  molecular  weight  of  the  light  gas.  The  simplified 
heuristic  mathematical  analysis  of  the  physical  phenomena  occuring  in  the  system  is  discussed 
in  the  following  paragraphs.  Table  1  defines  the  symbols  used  in  the  mathematic  analysis. 
Figure  2  illustrates  the  experimental  launcher  system  used,  and  Figure  3  shows  the  target  area 
and  target  evacuation  system. 

2.  COMBUSTION  TUBE  ANALYSIS 

The  solution  for  the  combustion  tube,  or  propellant  burning  side  of  the  device,  is  through 
a  standard  gun  ballistic  approach  similar  to  that  given  in  Reference  1. 

The  energy  balance  for  this  section  will  be 

E1  =  E2  +  e3  +  E4  (1) 

Where 


Ej  is  the  energy  put  into  the  system  by  combustion  of  the  solid  propellant. 

E->  is  the  translational  energy  of  piston. 

E3  is  the  heat  loss  to  waifs. 

£4  is  the  energy  required  to  accelerate  unburned  propellant  and  combustion  gases. 
The  chemical  energy  generated  will  be 

E1  =  mN  Cv  (TjJ  -  Tc)  (2) 

Reference: 

1.  Heiney,  O.  K„  Analytic  and  Experimental  Interior  Ballistics  of  Closed  Breech  Guns. 

Air  Force  Armament  Laboratory  AFATL-TR-69-42,  May  1969  (Unclassified). 


t  LAUNCH  TUBE 


Piston  wt:  152  grams 

Sabot  wt:  1  grain 

Powder  Charge:  300  •  BOO  grams 

Launch  Tube  inside  Diameter:  0.500  inch 


Pump  Tuba  Inside  Diameter:  2  300  inches  Launch  Tube  Effective  Length:  -132.5  inches 

Combustion  Tube  inside  Diameter:  4.0  inches  Sabot  Shearing  Pi  ensure:  20.000  tb/m 

Pump  Tube  Effective  Length:  135.9  inches  Pump  Tube  Outside  Dimeter:  4  inches 

Combustion  Tube  Effective  Length:  -71.25  inches  Launch  Tube  Outside  Diameter:  inches 


Figure  1.  Light  Gas  System 


TABLE  1.  LIST  OF  SYMBOLS  USED  FOR  MATHEMATICAL  ANALYSIS 


Bore  area  of  combustion  chamber 

Acceleration  of  piston 

Bore  area  launcti  tube 

Acceleration  at  payload 

Specific  heat  of  gas  at  constant  volume 

Total  mass  of  gun  propellant 

Force 

Acceleration  value  of  gravity 

Impetus  of  propellant 

Pseudo-mass  of  compression  piston 

Mass  of  compression  piston 

Molecular  weight  of  light  gas 

Mass  of  propellant  burned 

Mass  of  payload 

Mass  of  pump  tube  gas  charge 

Pressure  in  combustion  chamber 

Pressure  in  pump  tube  chamber 

Net  action  pressure  on  compression  piston 

Pressure  at  projectile  base 

Universal  gas  constant 

Linear  propellant  burning  rate 

Burning  surface  of  propellant 

Time 


10 


TABLE  1.  CONCLUDED 


Gun  propellant  gas  temperature1 
Gas  temperature  in  pump  tube 
Isochoric  flame  temperature  of  propellant 
Volume  of  pump  tube  chamber 
Velocity  of  compression  piston 
Initial  volume  of  propellant  chamber 
Payload  velocity 

Reference  distance  to  compression  piston 

Payload  distance  reference 

Heat  loss  factor  combustion  chamber 

Heat  loss  factor  pump  tube 

Gun  propellant  gas  density  factor 

Density  of  propellant 

Specific  heat  ratio  of  combustion  gases 

Specif  heat  ratio  of  light  gas 
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The  translational  energy  of  payload  will  be 


e2  = 


1 

2 


(3) 


The  heat  loss  of  the  gases  is  proportional  to  the  distance  traveled,  which  is  roughly  pro¬ 
portional  to  the  square  of  the  velocity  (Reference  2).  This  heat  loss  can  then  be  approximated 


e3  *  -Wc2 

4 


(4) 


Using  a  Kent  form  solution  (Reference  3)  with  high  velocity  modifications  for  the  energy 
contained  in  the  accelerating  gases  and  unburncd  propellant,  the  following  approximation  can 
be  obtained: 


E4  = 


(51 


In  this  equation,  6  equals  3  at  low  velocities  but  increases  at  high  velocities  because  the  density 
distribution  becomes  less  uniform.  This  effect,  and  the  variation  of  >'  with  payload  velocity, 
is  discussed  in  Reference  4. 


An  effective  mass  may  be  defined  as 


ma  =  mB 


+ 


CW 

9;; 


(6) 


Then 

E2  +  E3  +  E4  =  (1  +  pc)  -  J-  mavc2  +  PCACXC 
The  term  y  is  defined  by 


(vc-  1)  = 


Fe 

cvTo 


(7) 


(8) 


References: 

2.  Hirschfelder,  Kershner,  and  Curtiss,  Interior  Ballistics.  Volumes  I  and  II.  NDRC  Reports 
A-142  and  A-180,  February  and  April  1943,  (Declassified). 

3.  Kent,  R.  H.:  "Some  Special  Solutions  for  the  Motion  of  the  Powder  Gas,"  Physics  7,  1936. 

4.  Heiney,  O.  K.:  "A  New  Computer-Oriented  Formalism  for  Gun  Ballistics."  Proceedings 
3rd  ICRPG-AIAA  Solid  Propulsion  Conference.  Volume  I,  3-5  June  1968  (Confidential). 
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Then,  from  Equations  (2),  (7),  and  (8) 


mNpB  (-L— Ifi-)  =  -J-  (yc  -  1)  (1  4  /3c)mavc2  (9) 

The  temperature  ratio  is  eliminated  by  the  introduction  of  the  equation  of  state  to  give 
the  basic  ballistic  equation  for  the  propellant  combustion  chamber. 

W>  +  AcXc>  =  mNFB  -<Yc-  1)  (1  +^c>  ^2~  vc2  HO) 

The  following  differentia!  form  of  Equation  (10)  is  more  convenient  for  incremental 
computation: 


dPc 

dt 


(Vp  +  ACXC) 


dmN 

dt 


-  1)  (1  +  /fc)ma 


dt  dt 


p.  A  ^ 


LMc 


dt 


(ID 


Equation  (11),  coupled  with  the  expressions  for  propellant  burning  rate  and  those  describing 
the  motion  of  the  pro  ectile,  provides  a  complete  solution  for  the  combustion  chamber.  The 
expression  for  gas  gent  ration  is  then 


dmN 

dt 


rSg  P  p 


(12) 


where  Sg  is  the  total  exposed  propellant  burning  surface,  P  p  is  the  density  of  the  propellant, 
and  r  is  the  linear  burning  rate  of  the  propellant.  This  burning  rate  is  a  non-linear  function 
of  the  combustion  chamber  pressure;  thus,  the  r  vs.  P  data  must  be  read  into  the  computer  in 
tabular  form. 

The  equation  of  motion  for  the  compression  piston  is  derived  from  a  simple  force 
balance: 


F  -  pNac  “  mBap  *13* 

dvc  .  ,  _  VPc-pL> 

dt  p  mB 
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The  preceding  discussion  provides  the  solution  for  the  pressure  on  the  combustion  side, 
while  the  following  discussion  develops  the  solution  for  the  pump  tube.  The  solutions  are 
quasi  independent  and  coupled  only  through  piston  motion. 

3.  PUMP  TUBE  ANALYSIS 

The  solution  starts  with  an  equation  of  state  for  this  chamber; 


PLVA  =  nLRTL 
Then  taking 

R  =  Cv(rL  -  1) 

and  differentiating  gives 


(14) 


(15) 


dpL  _^a_  +  _  !jl_ 

dt  ■  1  '/  L  -  1  dt 


nLCV 


(16) 


Equation  (16)  is  a  differential  energy  equation  for  this  system.  Collecting  terms  and  including 
a  term  for  the  work  performed  by  payload  projectile  acceleration  gives 


dP,  R  dT, 

dt  L  VA  dt 


dVA 

dt 


-!>(_-  1)  ms  vg 


dVg 

dt 


(17) 


Consider  the  last  terr  with  mg  (handled  as  for  the  combustion  on  tube  analysis)  as  a 
combination  of  the  launched  payload  sabot  and  projectile  weight  plus  a  varying  fraction  of  the 
compressed  gas  mass.  The  advantage  of  the  use  of  helium  for  the  working  fluid  is  apparent 
here.  Figure  4  gives  the  value  of  this  variable  (gas  density  gradient  factor  0)  as  a  function  of 
projectile  velocities,  with  nj_  being  the  charge  of  helium  on  the  light  gas  side. 

m£=m  TJlk_  (18) 

»  g 


After  including  ^2,  the  result  for  this  term  is 


(rL  -  i)  (i  +  ^2)  msvs 


dt 


(19) 


This  gives,  finally,  the  basic  differential  pressure  equation: 


14 


■  O'  l_  -  1)(1  +P2)  msvs 


dtL  +  PlA: 
dt  VA 


dt 


dt 


dvS 

dt 

(20) 


The  result  of  Equation  (2)  is  then  coupled  with  the  results  from  the  combustion  chamber 
analysis.  Also,  required  is  the  following  relation  from  Reference  5  which  is  valid  to  the  first 
order  for  the  stepwise  quasiisentropic  approach  used; 

:  _P_  _riT_  =  _R _ dP  (21) 

T  dt  Cp  dt 


Again  using  Equation  (15)  gives 


(22) 


Equations  (20)  and  (22)  provide  an  incremental  solution  for  the  pump  tube  pressure  as 
a  function  of  time,  when  coupled  with  the  solution  for  compression  piston  motion. 

4.  PRESSURE  GRADIENT  AND  PAYLOAD  MOTION 

Equation  (2)  provides  a  time  history  of  the  space  mean  static  pressure.  An  expression 
for  payload  motion,  however,  requires  the  pressure  at  the  base  of  the  shot  to  be  defined. 
Reference  1  covers  this  pressure  gradient  computation  in  some  detail.  The  results  of  that 
analysis  provide 


(23) 


Where  Pg  is  pressure  at  the  projectile  and  P|_  is  mean  chamber  pressure. 

Equation  (23)  is  plotted  in  Figure  5  and  dramatically  illustrates  the  advantage  of  using  a 
low  molecular  weight  gas  as  the  driving  medium.  It  is  seen  that  at  a  velocity  of  10,000  fps, 
combustion  gases  with  a  molecular  weight  of  24  are  no  longer  able  to  deliver  energy  to  the 
accelerating  payload.  Helium,  however,  with  a  molecular  weight  of  4  is  seen  to  remain  44 
Reference; 

■5.  Liepmann,  H.  W.,  and  Roshko,  A.:  Elements  of  Gas  Dynamics.  Wiley,  1957. 
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percent  efficient  in  it.  ability  to  deliver  useful  energy  to  the  ^projectile.  At  these  velocities, 
regardless  of  the  breech  pressure,  a  gas  with  a  he^-vy  mean  molecular  weight  of  24  will  require 
virtually  all  the  generated  thermal  eneryy  to  acquire  the  necessary  directed  kinetic  energy 
to  reach  payload  velocity.  The  lighter  molecular  weight  gas  will  require  much  less  of  the 
available  energy  tea  be  accelerated  to  payload  velocity.  With  this  pressure  gradient  defined, 
the  payload  motior  equations  become,  then,  quite  simple  in  an  incremental  form: 

vS  =  VS  +  aS  * 1  i  (24) 

t 

I 

with  ! 

) 

as  "  pSAS'/mS  j  I25* 

f 

I 

Equations  (11),  (20),  and  (24)  provide  solutions  for  pressure  in  the  combustion  tube,  pump 
tube,  and  payload  motion,  respectively.  A  complete  solution  is  possible  if  the  following  initial 
conditions  are  known:  volumes,  diameters,  and  weights  of  projectiles  and  propellants.  Addi¬ 
tionally,  initial  gas  pressurization  and  sabot  release  pressures  must  be  determined  accurately 
because  the  system  performance  is  extraordinarily  sensitive  to  these  two  parameters. 


SECTION  III 


USERS  GUIDE 

This  section  provides  a  nomenclature  table  (Table  2),  a  program  listing,  and  sample  input 
and  output  data  so  that  this  computer  code  can  be  run  as  it  currently  exists  or  it  can  be 
modified  for  other  particular  applications,  such  as  those  discussed  in  Reference  6. 

An  abbreviated  flow  chart  is  provided  in  Figure  6  and  provides  general  comment  statements 
that  may  be  lacking  in  the  program  itself.  The  program  may  be  run  on  reasonably  small 
computers  because  core  requirements  are  moderate.  Table  3  provides  a  sample  of  the  necessary 
case  input  data,  and  Table  4  is  the  associated  output 

1.  INPUT  DATA 

The  input  data  for  the  program  are  in  two  main  categories.  The  first  goes  in  only  once, 
and  consists  of  13  cards  containing  propellant  and  light  gas  data.  The  second  set  of  data 
consists  of  case  cards;  three  cards  are  required  per  case,  arid  as  many  case  cards  may  be 
stacked  as  is  desired.  The  particular  data  are  as  follows: 

a.  Propellant  and  light  gas  data: 

Card  1  contains  propellant  impetus,  specific  heat  ratio,  density,  covolume,  and  type. 

Cards  2  and  3  contain  20  reference  pressures  ftr  propellant  linear  burning  rates. 

Cards  4  and  5  contain  20  burning  rates  at  the  fixed  reference  pressures. 

Cards  6  and  7  contain  20  fixed  propellant  gas  velocities. 

Cards  8  and  9  contain  20  propellant  gas  density  factors  corresponding  to  mean 
density  distribution  in  systems  with  the  fixed  propellant  g?s  velocities. 

Cards  10  and  11  contain  20  fixed  helium  gas  velocities. 

Cards  12  and  13  contain  20  system  helium  density  distribution  factors  corresponding 
to  the  fixed  helium  gas  velocities. 

b.  Case  Cards: 

Card  1  consists  of  the  following  system  physical  property  data  for  the  combustion 
chamber.  Bore  area,  chamber  volume,  piston  weight,  piston  travel,  propellant  web,  heat 
loss  factor,  propellant  charge,  initial  helium  pressure,  piston  shot  start  pressure. 


Reference: 

6.  Rynearsori,  R.  J.:  Optimization  of  a  Two-Stage  tight  Gas  Gun.  Thesis,  Texas  A&M 
University,  December  1972. 
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TABLE  2,  PROGRAM  NOMENCLATURE 


| 


S  •  ACEL 

Acceleration  of  the  piston 

i 

i 

i 

AvLEA 

Cross-sectional  area  of  barrel 

J 

ALGT 

Acceleration  of  the  projectile 

] 

>  t 

AREA 

Cross-sectional  area  of  pump  tube 

j 

!  BARE 

Burning  area  of  propellant 

J 

j 

j  8DIS 

Piston  travel 

■j 

1  BETA 

Heat  loss  coefficient  for  propellant  gas 

BL!S 

Projectile  travel 

BLTA 

Heat  loss  coefficient  for  light  gas 

BP 

Mass  of  propellant  burned 

CHG 

Initial  powder  charge  weight 

\ 

i 

!  CPEG 

Projectile  base  pressure 

j 

I 

1 

i  CPIC 

1 

Projectile  muzzle  velocity 

i 

i 

j 

! 

j  CPLG 

Average  light  gas  pressure 

j 

i 

!  CPRS 

j 

Average  propellant  gas  pressure 

\ 

,1 

1 

1  CVL 

Oavolume  (’ll  for  light  gas 

< 

j 

i 

j  CVOL 

Powder  chamber  volume 

\ 

] 

;  DAN 

j 

Diameter  for  ball  type  propellant 

1 

]  DELTA 

Time  increment 

I 

1  ,  DFLG 

",3"  factor  for  light  gas 

„  j 

i 

DIN 

Propellant  grain  inside  diameter 

j 

j  '  DLDT 

Piston  velocity 

i 

]  DOT 

Propellant  grain  outside  diameter 

i 

i 

j  DPDT 

Pressure  slope  in  propellant  gas 

i 

I 


1 


I 

i 


I 


DPLG 

DTDT 

EFM 

FIMP 

FPU 

F  VOL 

GAMA 

GIN 

GMLG 

GN 

GRNS 

HCPP 

HGBL 

HGID 

HGIP 

HGIT 

HGIV 

HGM 

HG3M 

HLGP 

IPT 

PIT 

PFAC 

PLGR 


TABLE  2.  CONTINUED 
Pressure  slope  in  light  gas 
Temperature  slope  in  light  gas 
Effective  or  psuedo  mass  of  piston 
Impetus  of  propellant 
Fraction  of  propellant  burned 
"Free"  volume  of  propellant  gas 
Specific  heat  ratio  of  propellant  combustion  gas 
Initial  powder  gas  density 
Specific  heat  ratio  of  light  gas 
Propeliant  gas  density 
Number  of  propellant  grains 
Peak  pressure  in  propellant  gas 
Barrel  length 
Initial  light  gas  density 
Initial  light  gas  pressure 
Initial  light  gas  temperature 
Initial  light  gas  volume 
Mass  of  light  gas 
Weight  of  projectile 
Light  gas  peak  pressure 
Identifies  propellant  type 
1.5,  empirical  correction  factor 
"6"  factor  for  propellant  gas 
I  sen  tropic  pressure  ratio  in  light  gas 
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TABLE  2.  CONCLUDED 


PREX  Piston  base  pressure 

PRS  Array  of  powder  chamber  pressures 

PSL  «'•  array  for  tight  gas 

PSY  .V  array  for  propellant  gas 

RFST  Diaphragm  burst  pressure  or  projectile  release  pressure 

RHO  Propellant  weight  density 

RTF  Universal  gas  constant  (Units  =  ft  •  lbf/lbm0|e°K)  value  =  2780 

RUN  Pump  tube  length 

SABPR  Powder  diaphragm  burst  pressure  or  piston  release  pressure 

SCPRS  Initial  powder  chamber  pressure 

SHOT  Piston  mass 

TIME  Elapsed  time 

TF  Flame  temperature  of  propellant 

TLGS  Average  light  gas  temperature 

TYPE  Output  array  to  write  propellant  name 

UBW  Unburned  propellant  volume 

VEE  Velocity  array  for  propellant  gas 

VEL1  Old  piston  velocity 

VEL2  New  piston  velocity 

VLE  Velocity  array  for  light  g3s 

VLG  Projectile  velocity 

VLGS  List  gas  volume 

WEB  Propellant  grain  thickness 

WMAL  Light  gas  molecular  weight 

WMOL  Propellant  gas  molecular  weight 

XLIN  Length  of  propellant  grain 
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TABLE  3.  INPUT  DATA 


346180, 

1.232 

.0603 

29,69 

M-10 

300. 

300. 

700, 

1000. 

1500. 

2000. 

2500. 

3000. 

4000.  5000 . 

6000, 

8080 , 

10000, 

20000, 

30000 . 

40000. 

50000. 

70000. 

100000,200000, 

.13 

.2 

.28 

.38 

.32 

.68 

.81 

.96 

1.2  1.43 

1.72 

2.20 

2.70 

4.6 

6.7 

8.3 

9.6 

12.2 

15.8  25. 

0. 

500. 

1000. 

1500. 

2000. 

2500. 

3000, 

3500. 

4030,  4500. 

5000 , 

55fe0, 

6000. 

6300, 

7000. 

7  500, 

8000, 

8500 , 

9030.  12000. 

3. 

3. 

3.05 

3.1 

3.  IB 

3.35 

3,38 

3.50 

3.70  3.85 

4,03 

4.3 

4,55 

4,65 

5.2 

5.55 

5.95 

6.40 

7.03  12.0 

0. 

4000. 

’6000". 

8000. 

10000, 

V2000. 

1 400P , 

16003.  18020. 

20000. 

22000.  24000. 

26000. 

28000  . 

32000. 

36000. 

40000 . 

45003.  50000. 

3. 

3. 

3.  1 

3.2 

3.4 

3.6 

3.9 

4.2 

4.6  5.1 

5.3 

~  6.1 

6.7 

7.3 

8.0 

9.6 

1  1.5 

13.6 

16.7  20.2 

4,15 

895.  .45 

136 

.  .018  » 3 

.32 

100. 

900. 

3000 

.  132. 

40. 

0033 

565. 

300. 

• 

196  .2 

5 mCT  «T 


£-«9Cf 


EB 


TABLE  4.  CASE  OUTPUT 

C*B  B3«E  aPEa 


. .  *.45 

?.52P  7, 

?182  135. 

7 

""8  5  3.  pa 

. . J.t3 

LI3-*  r,*S  G 

'JV  c  4  T  * 

9-ct  -etoT 

B .LtN&T* 

c *e  vcte. 

ec«t  area 

"til  tISS 

— . 

?.??33e 

132.  ?  ? 

365.22 

2.196 

?,?? 

•‘CuECOL4»  mf 

rS»T  •  21.? 

-E4T  LCSS 

f  * 

rw  JS 

2.32  p«cpe 

;,SS3 

s»#rt“  is  ».ia 

POOPEWL**'’ 

fte*  !S  S!sG 

pcatc^m 

r  * 

_CC*5? AST 

3US?  *C£ 

_ _  ..  _ 

_ 

•••«••••*•« 

Bi*DP£*wi.4S?  S3 

?! 

*  •  .  .  . 

G-T  545  4IDE 

TS“t 

C*a»*9  saES 

T3iVEu 

»92P  Bt'R^E 

:*«E5  si:*e 

ve_:c:t« 

•*:s  PRES  >EvC£ITv 

*e*»Ew 

93  *5E5 

*f  ■» 

P.70P 8??? 

is?. 5 

a.ea* 

*,2116 

21383.27 

2,?? 

4?.??  P.00 

0,0? 

4?.< 

3??.?2 

P.7216??? 

13*,! 

?.a?? 

2,2231 

?i  384,44 

?.2P 

4?.??  ? , ?  a 

2,2? 

4?,  ?? 

”579224??? 

IBS  ,S 

?.??? 

2.0347 

’  2 ! 325, 62  ■“ 

P.PP 

*?.»  >.n 

P.PP 

-<?.?? 

P.P23277P 

its.  i 

a. 2/p 

2.2482 

21386.79 

2.2? 

4?.??  7.22 

P.7? 

47,02 

3P7.72 

?,??4??2? 

18i,3 

?,??? 

2. 0376 

21587,87 

?./? 

4?.r?  ?,2? 

0.2? 

42  ,?7 

322 .22 

p, 20402?? 

2?3 , 5 

a.??? 

2,0883 

21389,1- 

/,?2 

4?,2?  ?,00 

2.2? 

42.2? 

370,?2 

P,PP567?? 

22  ? ,  9 

?.??? 

2,2908 

21582,32 

e.ea 

4?,??  7.0P 

2.?? 

42.02 

320. 77 

P.7764P2? 

23!.? 

a.??? 

2,0324 

21391.3? 

?.?? 

4?,??  P.0? 

e.?a 

42,77 

32?, ?? 

P,*8Jiaa:> 

255. 4 

j.j?* 

!. iei« 

'  ? 1 3 9i , Sfi 

?.P7 

•if, it  l.fi 

p.ii 

4?  ,?2 

37  0  >2 

p.pae???? 

2?2.7 

a.??? 

0.1136 

21393.83 

i.n 

42.02  ?.?? 

e,7? 

47, ?P 

3? 2 , 2  ? 

2.70860?? 

2  5  S  ,  P 

?.??? 

2,1271 

2 1 333,03 

2,2? 

4?,??  ?.?? 

e.2? 

47.2? 

3??,?? 

p,*096??? 

3?7 , 3 

?.  ??? 

3,1367 

21893,62 

7. It 

4?,??  7.22 

P.P? 

42.7P 

3? 7,2 2 

P.01P422? 

323.3 

?.??? 

3,1327 

22934.52 

?.?? 

42.??  0.22 

P.72 

47.22 

322,72 

P.PU2??? 

3*4,1 

a.??? 

3.1633 

24224.78 

7.2? 

4?.??  ?,?? 

e.p? 

4  2.2? 

?;  pi????? 

383.2 

2,1*53 

23167.04 

?.?? 

i?,??  a,?? 

?.?> 

47,7? 

3P7  ,?? 

Jii.i 

i  ,z?? 

3,1923 

76353,73 

7,02 

4?  .  ??  ?,2P 

P.2? 

47,22 

370^22 

4  ?  ft  .  ? 

?,??? 

2,224* 

27917,34 

?.2? 

47.??  7.2? 

P.72 

42. PP 

3?2 , 7? 

P  .5  i  44???. 

428.5 

?.??? 

a,?2?2 

28931,*  5 

?  .  ?  3 

42.2?  ?.?? 

P.77 

At.  ?P 

2?r J?? 

t. until 

4  22,! 

?.??? 

3.2339 

30327,44 

7, 2? 

4?  ,  ??  ?,27 

2.72 

42.2? 

*.»!!»!.’» 

4  7  7.5 

?.?,;? 

3,2523 

31749,45 

?  ,2? 

42.0?  P.27 

p,a? 

322, 7? 

?.?( 58*?? 

57  3.3 

2.??? 

0.7598 

■55262.35  ' 

2.3? 

:i .??  ?,37 

2.;p 

4?. 22 

!??’» 

p .7 1 ?9»p? 

331.? 

2.77? 

2,298,3 

35?75. 34 

3.2? 

4?.??  ?,?? 

?.?? 

42.2? 

3>?’?2 

? 64??.? 

259.8 

».  ??? 

2,307* 

36989.99 

2.?? 

4?.??  ?.?? 

2,02 

42.72 

P.P192??? 

•5?,  4 

?.??? 

2.326C 

55008,2? 

p.e? 

4?.??  0,30 

e.?2 

47,?? 

3?P^?2 

7 ,727777? 

522.3 

?.??? 

2,3453 

41136.96 

?.?? 

47. ?2  P.2P 

2.?? 

42.2? 

?.???«??? 

a.??? 

3,3721 

43382,36 

t.tr 

47.7?  e.2? 

2.2  2 

7,021607? 

fit. 2 

?.??? 

'P.3360 

4375?. B7 

l.fr 

..as 

.  "2.2? 

42.2? 

7.722477? 

7  2  9,9 

?.*/0 

2.4211 

47967.18 

?.?? 

*?.>?  r.tt 

2,22 

4?  .2? 

P.7234??? 

769.2 

?.??? 

2.4474 

inti.n 

a. 

4?.??  ?,20 

?,PP 

7,?2*??0? 

Si*. 3 

a.?2? 

3.4748 

32418.22 

7.?2 

<?.>.•  p.aa 

0.22 

42.0? 

7. 9243??? 

833.3 

a.??? 

3.3236 

34797.26 

2.0? 

<?,??  r.tt 

?  .  22 

42.2? 

7,7236227 

886,3 

a.??? 

3.3336 

37294. ge 

0.2? 

4?,ee  ?.?? 

2.72 

?.?s))|2:c 

5?8,1 

?.?22  ** 

a;34?5 

’3?8it,l5 

la.sr 

it.il  i>.*i  • 

"  2.0? 

42.  ?2 

7,72377-7 

9?fl  .  4 

a.ai? 

e, 344? 

37980, !3 

2?, 36 

4  ? » ?  s  ?.a? 

2.2? 

*  4?,  21 

3?  2  ’.*2 

7.72S61S7 

8l?.7 

?.fl22 

3,3456 

37399.41 

32,89 

40.es  p.4i 

2.02 

».f S385JJ 

«IS.J 

i.iK 

3.5471 

37438.74 

41,23 

40, ?3  e.P0 

a  .a? 

7. *25855? 

3  !  2 , 3 

P.032 

3.5467 

3?307, *3 

31.6? 

<?,n  r.n 

2.70 

41,** 

7,72383*7 

917,5 

2.798 

J 

57i?2.83 

62, ?? 

4?,  06  ?.?? 

P.2? 

*,a250?3. 

515.5 

7.121 

8.3318  ' 

'37034. 0»  -  ” 

"72.13 

il.il  ~  t,?3 

"  “P.PP 

4P.fi 

SPP.23 

*,*j4*fS* 

*22,1“ 

•  .IS* 

924.4 

0.201 

0,0240950 

929,7 

0.248 

«*.  1 33i* 

929,9 

0.331 

0, *281730 

931,2 

0.358 

0. 028213? 

633,4 

0.421 

0,0262530 

935.7 

*.488 

0.0282950 

637.9 

2,561 

0.1*263350 

940.1 

0,839 

SK 

942,4 

8.722 

P.02841S0 

544,  * 

0,810 

0.026455? 

646.6 

0.9?3 

0.0284950 

649,? 

1,001 

0.028333* 

931,2 

1.103 

0, *28373* 

933,4 

1.213 

63  3. 3 

1.32? 

0.02665'J 

957.7 

1,445 

*.0288949 

959,9 

1.371 

0,0297349 

962,0 

:.7y0 

0.0287740 

9  6  4,? 

1.833 

0. 0*68149 

966,3 

1.973 

B, **88344 

698.4 

2,120 

0. 0286949 

670,6 

2.271 

0.0269349 

07?,? 

2.427 

0,0299749 

974,8 

2.  588 

0,0277149 

576.9 

2.734 

0. *270549 

979 ,0 

2.926 

0,0270949 

981,0 

3,103 

O.O20I349 

9b  3,  l 

3.285 

0.0271749 

085.2 

3.473 

0. *272149 

687.9 

3,666 

0,0272549 

9  8  5.3 

3*864 

0.02 72949 

661.3 

4,088 

"0.027334? 

9  9  3.3 

4.277 

0.0273749 

963.J 

4,492 

0**274149 

997.3 

4.711 

*.*274349 

969,3 

4,937 

O.OJ74949 

1001,3 

3,187 

0.0273349 

m  3,3 

3.403 

8,0275749 

1275, 2 

5,  645 - 

0.0276149 

too? ,? 

3,882 

0.0273549 

1 

6,144 

0.0276949 

1*11,0 

6.402” 

0.0277349 

1012,9 

8.665 

0.0277749 

1014,8 

8. <34 

0.0276149 

1215.7 

7,278 - 

0,0276349 

1019,6 

7,468 

0.0276649 

1020.3 

7,773 

0,0279349 

1022.3 

8.064 

0. *>279749 

1024.2 

6.380 

*.»28f  !*■> 

1024.0 

8,682 

~*,*2!j349 

1*8?.! 

8,969 

0.0280949 

1*29,6 

9.281 

0.0281349 

1031.4 

9,600 

0.0291749 

1033, 2 

9.923 

0.0282149 

1033.0 

10.233 

0,0262349 

1036.8 

10.307 

T. *282949 

1038,5 

~!r,'S2f — — 

0.0283349 

1040.2 

1  1.274 

TABLE  4.  CONTINUED 


0,3534' 

56892.96' ' 

62,88 

40,  U 

0.5550 

56747,80 

93.36 

40.  M 

0,3566 

56599,17 

183.87 

4*.  17 

0.6381 

36447.10 

1H.48 

M.fi 

0.3597 

56291.39 

124.96 

<0.25 

0.3613 

36132.68 

133,54 

*0.30 

0.3529 

33970.36 

146.13 

4fe,35 

0,5645 

33804,67 

136,78 

48.40 

0.3661 

33635.63 

167..J3 

40,46 

0,5677 

13463,25 

176.12 

40-;s2 

0.5693 

55287.55 

186.82 

40.36 

0,3709 

55108,56 

i 99,53 

40,63 

0.3725 

34926.31 

212.3? 

4*. 72 

0.3741 

34740.80 

221.29 

40.60 

0,5758 

34552,09 

231.68 

40.88 

0.3774 

34380.10 

242.78 

40; 96“ . 

0.5790 

54165. 11 

233.33 

41,05 

0,5806 

53966,90 

264.42 

41.13 

0.8*23 

53765.59 

273.31 

41.24 

0.5839 

53361.21 

285.23 

41.34 

0,5853 

53353.78 

297.16 

<1,43 

0.5872 

53143.35 

328.12 

ir;** 

0.5888 

32929.94 

319,10 

41.67 

0.5903 

32713.38 

33?. 10 

‘1.79 

8.3921 

52494,32 

341.12 

41,92 

0.3938 

32272.18 

332.16 

4?, 04 

0.5954 

•  52047.22 

363.23 

42.18 

0,5971 

51019.45 

374.31 

‘2, 3i  — • 

0.5987 

51368. 9i 

365.41 

*2.43 

0,6004 

51355.69 

396,53 

<2,5e 

0.8021 

51119.76 

407.68 

4?.  73 

0  ,6036 

50681  ,n 

*16,94 

*2.91 

0.6034 

50640,05 

43.*. *2 

43,07 

0.8071 

50395.34 

4*1.2? 

"  '  ‘"43723 - 

0.6088 

5*158. 13 

432.43 

‘3.41 

0,6105 

49901 , 45 

463.67 

43,38 

8 .6 1 22 

49650,34 

474,6? 

*3.76 

0.6138 

49387,04 

486.19 

43,93 

0.6153 

49114.94 

497.48 

44,14 

0,6172 

40848. 2B 

508.7* 

427  3  4 

0,6169 

48363.88 

520.11 

44,34 

E  ,620  6 

48285.45 

531.44 

44.73 

0.6223 

48005.02 

542.8? 

4.1.96 

0.6240 

47722.67 

634.1? 

43,16 

0.6237 

47438.43 

365.33 

45,41 

3.6274 

. 47152.35  • 

‘  575,93 

45-min - 

0.6292 

46864,51 

388.35 

45.88 

0.5309 

46574,94 

699,79 

46,12 

0.6326 

46283.71 

611.24 

46.37 

0.6343 

45990,86 

622.69 

46,63 

0,6350 

455P6.46 

634,16 

<6.89 

0,  K377- 

43400.55 

645, 55 

■  “  4?;T5" 

0.6395 

43103,20 

657.15 

47.44 

*.6412 

44804.46 

668.65 

47,72 

0.6429 

44584.39 

6B0.18 

48.  *1 

0.6447 

44203.04 

891.71 

48,31 

0,6464 

43900,48 

703,23 

48,61 

4339^,76 

rrr.fl ! 

ISTW 

0.6499 

43291.94 

726.38 

49,24 

0.C0 

0.ea 

40.lt 

300.33 

0.00 

6.00 

40.14 

300,42 

0,00 

0,00 

4  0  a  1  7 

300,32 

0,00 

0.00 

40,21 

300.63 

0.00 

e.00 

40,23 

300,76 

0.00 

f.Biff 

*e.30 

300,89 

0.00 

0,00 

40,33 

301.04 

0.00 

0,00 
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C4JE  IfRMIWiTEn  tS  MA*IHU"  VELOCITY  -REACHED-  AT  '  35.4  IWCHE5  OF  TRAVEL 
LIGHT  GAS  PEA«  PRESSURE  15  17750.9  P$I 
COHUSTIOR  SIDE  PEa«  PRESSURE  IS  1J46.4  P5I 
LIGHT  G4S  Guh  “uIIlE  vElOCITt  IS  I l 360 . 8F T/SEC 


Card  2  consists  of  pump  tube  physical  parameters,  which  are,  in  order:  sabot  start 
pressure,  launch  barrel  length,  initial  pump  tube  pressure,  payload  weight,  pump  tube  volume, 
pump  tube  initial  temperature,  launch  tube  area,  and  pump  tube  heat  loss. 

Card  3  has  a  propellant  form  flag;  the  number  1  indicates  a  constant  burning  surface. 

2.  PROGRAM  OUTPUT 

The  program  listing  is  given  in  Table  5.  The  output  (Table  4)  is  for  the  listed  input  data 
(Table  3)  and  corresponds  to  an  experimental  test  firing  conducted  with  a  light  gas  gun.  The 
correlation  of  experiment  to  analysis  for  both  the  pump  tube  and  the  combustion  tube  are 
given  in  Figures  7  and  8.  The  time  sequencing  for  the  analytic  data  has  been  worked  backward 
from  the  experimental  peak  pressure  because,  in  practice,  the  ignition  transient  is  exceptionally 
long.  The  phenomenally  long  ignition  delay  is  due  to  the  extremely  low  propellant  loading 
density  and  attendant  required  flame  spread  time  at  low  pressure.  Because  of  the  erratic 
combustion  properties  of  gun  propellant  at  very  low  pressure,  this  several  hundred  millisecond 
ignition  delay  is  effectively  non-analytic. 

The  computer  output  shows  that  the  program  will  print  data  every  800  microseconds  until 
the  driving  piston  starts  to  move,  in  this  case  at  900  psi.  The  program  then  prints  data  at 
a  rate  of  every  40  microseconds  until  the  payload  sabot  is  sheared  and  starts  to  move.  This 
occurs  at  5,000  psi  light  gas  pressure,  and  the  final  print  frequency  rate  becomes  once  every 
10  microseconds. 

Figures  7  and  8  present  a  correlation  of  experiment  and  analysis  for  the  comt  ustion  side 
and  for  the  pump  side,  respectively.  As  previously  mentioned,  time  justification  ,ias  been  by 
coincident  pressure  peaking  on  the  pump  tube  side. 

As  is  seen,  the  correlation  for  both  the  helium  and  propellant  chambers  is  very  good; 
however,  it  is  not  exact.  The  pump  tube  side  has  a  somewhat  sharper  pressure  decay  experi¬ 
mentally  than  is  predicted  by  the  analysis,  this  is  doubtless  due  to  heat  loss  effects  or  sabot 
friction  which  is  not  totally  accounted  for  in  the  mathematical  model.  The  combustion 
tube  pressure  rise  has  an  analytical  variance  from  experimental,  which  is  to  be  anticipated  from 
the  non-analytic  nature  of  the  slow  smouldering  ignition  transient  described  previously. 
Additionally,  at  the  end  of  the  compression  stroke,  the  analysis  predicts  a  slight  pressure 
recovery  ,  which  has  not  been  seen  experimentally.  The  analytical  muzzle  velocity  prediction 
of  11,360  ft/sec  compares  with  ’he  experimental  velocity  of  11,210  ft/sec. 

In  summation,  it  can  be  stated  that  apparently  the  mathematical  mode!  of  the  two-stage 
light  gas  system  is  valid  for  performance  prediction.  It  serves  as  a  useful  tool  in  determining 
required  propellant  and  gas  charge  loads  and  in  optimizing  sabot  and  piston  design  shot  start 
pressures.  Although  useful  refinements  should  be  readily  apparent,  the  core  of  the  program 
can  be  easily  manipulated  to  accept  more  sophisticated  modeling,  if  required.  As  a  heuristic 
design  tool,  the  core  program  is  largely  satisfactory. 
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TABLE  5.  PROGRAM  LISTING 

DIMENSION  P9Y(20),VEE(20) 

01  MENS  1 QN  PR  9(2?} ,B*TE  (2*) ,!YPE (8) 

01  PENSION  prop?  (4,31 

—  dimension  rsrn  s  i  vtmt rroTTHTZTi tj  tottatt) - 

DIMENSION  P9L  (20} , VUk (20} 

COMMON/  /PCPRS U 000 ),PBO IS (1000}, RUN, P?OP, NPTS 

"  COHHON/PLTHDR/  SHOT, CHG, NEB, CVOL, AREA, VH  . - 

RE*0(4,3}FlMp,0*Mi,RhO,CVL, (TVPE(N} ,N*l,4J 

3  F0RM*T(FP.l,F8,3,f 8.4,Fe.2, 4 42) 

- "rfadtavoprs-  "  - - ■- 

4  PORM*T(l0Fa, 1) 

RE*0(4,3)R*TE 

“  “  5  FORMAT  (1PF8 , 3)  ' 

RE  AO (4, 24) VEE 
24  FORM*! ( 1PF8, 1 ) 

- UMDN.JGJPSy - - - - 

2«  FORM*! (10F8.3) 

RE*D(4,24)VLE 

- READ  (4,2S5  PSL - - ““ . 

PIT*  1 , 5 

_ T*0*G*M*/(C*MA-1,) 

7  F0RH*T  (FT. 3, F7. 3, F7. 3,  F7.1.F7, 4^7, 2^7. 3, 2F8.I) 

RE *0(4, 14)  RFST,N0BU,H0IP,H59M,HGIV,HGIT,ALE*,8LT* 
“  l'<  FORH*T(BFl0,0) 

READ  (4, 110)  IPT,0*N,00T,0IN,!<LIN,MM01„ 

1  10  FORM*!  (11,6718,7)  _ _ 

16  NPTS »0  * 

EFM»0 

-  ACEI*0  . . . * 

VEL 1 «0 
AV0L4CV0L 


R*T40, 

NPTS«0 

KUV*S  '  '  . .  " 

Bt IS  40 , 

HLGP40 

"HCPPi*  . 

ILF*1 

JFG40 

RTP»27S0, 

KUE»-1 

VS.G40, 

^-augtw; - — - 

TUG34HGIT 

OTOT40, 

GMLGil .87 
SLG*GMLG-l. 

NN*L«4, 

"DPU5b0, - 

GAG*1.-(1./GMLG) 

CPLG4HGIP 
CPEG4CPUG 
T I ME»0  , 

OELTA*0.0001 

■VLG54HGrW— - - - 

HGID4MG!P*Mm*l/ (R!F*HGIT*12,) 
HGMiHGIO*HGIV 

BP»0. 

VEL40. 
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TABLE  5,  CONTINUED 

IF(WMQl.lE.0.0)  wnOL*24, 

TF»WMOL»FIMP/RTP 
X  TF«TF 

GIN»SCPRS«WM0L/(»TF*TF«12.) 

fvol*cvol«<chg/rho)  ~ 

GO  TO  (125,  123, 120, 123, 123) , IPT 
122  FPU»B. 

RT  «0  . 

SARI«3.«CHG/ (RH0*wE8) 

GO  TO  132 

~TJ3  AOToOOT.OOT. 3. 1*17/4.  . .  "  . . .  “ 

AIN»0IN«0IN*3, 1417/4, 

IF  (IPT.EO, 4)  AIN«7,0«AIN 
IF (IPT. £0,3)  *lM«19.0*4I«i 
*|F«*OT-*IN 

GRNSwCHG/(AEF«XLIN*RH0) 

- GO  TO  130 . . . . .  . . 

129  9AREn2,oCHG/(RHO*WEB) 

130  WRITE (8,8) 

8  FORmrtlHl.eiH  SHOT  wT,  CHARGE  WEB  B.LENGTH  CMB  VOL 

1BORE  AREA//) 

WRITE («, 9) SHOT, CHG, WEB, RUN, C VOL, AREA 

— 9  FORHAT[Fi0.2,Fi0.3,Fi2,4,narr,Fr0.T,n?,2/r/y - - - 

WRITE (6,2fl) 

29  FORHAT (24H  LIGHT  GAS  GUN  data  , 60H  SHOT  WEIGHT  B.LENGTH 
1 Ch8  VOl,  SORE  AREA  HEAT  LOSS) 

WRITE  (8, 25)  HGSM,HG8L,HGIV,ALtA,BLTA 
23  F0RHAT(27X,F11.3,F1  1.2, FU. 2, F  12.3, F 12. 2) 

-  OATAtU.BHOL  ’  ~  _  "  - - -  ' - - 

OATA(2)«XLIN 
DA  T  A ( 3 ) BOOT 
OaTA  (4 ) «OIN 

WRITE  C»,30) wmol, BETA, (TYPE (N) ,N« 1,4) 

30  FORHAT  C//21H  HOLECULAR  WEIGHT  •  ,F3.l,23H  HEAT  LOSS  FACTOR  T3  , 

- IF3.2.35H  SRCPELLFNT  USED  IN  "SYSTEM"  1ST  jnTJTT . "  . . 

GO  TO  (31, 33, 33, 37, 61), IPT 

31  WRITE(8,32) 

33  FORHAT (39H  PROPELLANT  FORM  IS  SINGLE  PERFORATE  OR  CONSTANT  SURF* 
ICE//) 

■  .'0  TO  80 

'"3T  *  irtET®,  34) . . . 

■>1  0PMAT(47H  PROPELLANT  FORM  13  DETERRED  SINGLE  PERFORATE//) 

GO  TO  60 

"33  WRITE  (8,36)  - - - ~ - - -  - 

38  FORHAT (33H  PROPELLANT  Form  IS  DETERRED  BALL//) 

GO  TO  60 

37  •••RITEC8.33)  -  - - - - - — - : - 

38  FQRhaT («7h  PRoPELIAn;  form  IS  OETERREO  SEVEN  PEHFQRATE  //) 

GO  TO  S3 

61  WRITE (S , 6J ) 

62  FORMAT (49H  PROPELLANT  FORM  IS  OETERREO  NINETEEN  PERFORATE//) 

60  WRITE (8,93) 

93  FORHir(2V,'Hff.WiWL*#**f«r,T0T,T8fl - PROPELLANT TTUr,  mVTTKWTrWT*  - 

. . . . .  GAS  S ! OE , 1 0X , 1 0H.««* «*•••• ) 

WRITE (f , 96) 

98  FOHMAT(8X,4MTIHE,4*,10HCHAMB  PRES, 4X, 6HTRAVEL , 4X, l 1 HPROP  BURNED, 1 0 
1HPRES  SLOPE, 3Y,  17HVELOCITY**CB  PR ES , 3X , BHvELOC I TY , 4 X , 6H TR A VEL , 3 X , 7 
2HB3  PRE3,6X,4HTEHP) 

- BP«0.0  '  - - 

vEL2*B« 

BDISW0.0 

PTOR»3000, 

CPRS»3CPRS 
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TABLE  5.  CONTINUED 

GAHBa(1.*8ETA)a(GaPA»1,) 

DO  39  J»  1 ,  3P*>0 ,  1 
IECBLIS.GE.hGBL)  GO  T0  39 

GO  TO  C74, 135, 74,134, 1341, IPT  .  _ 

— iTS  IFtRAT-AEBi  135,74,74 

135  RAT  nR*. 00085 
OINaOIRaRaT 
*LIN»XIIN-9AT 
AINlQlN.niN. 3, 1417/4, 

A1X«0I«»3,1417  _  _  _ 

~  IE~ClPt.Ni.4J  GO  to  1 32 
AINa7,0aAJN 
AIX«7,0«A1X 
GO  TO  138 

132  IE(IPT,NE.31GO  TO  136 

AIN»19.0aAIN  _ 

-  A  J  X  »  1"5“,  0  *  A  i  N 

136  AEEaaOT-AIN 

BAREaGRNSa {2,«A£EaXLIN**IX) 

""'74  J  A  ■  1 

IE(CPRS,LT»300«0)  GO  TO  76 

78  IF(CPRS-PfiStJA) ) 77,76,76  _  _ 

- TT  Ji.JT*!  “  ‘  . . 

GO  TO  78 
78  RaRATE  C J  A) 

GO  TO  10 

77  OITaPRS CJAJ-PRS [ jA-n 

OATARATEjyA).RATEtJ*-l)__  _ __ _ 

p  ig'bcprs-prs T  Ja-  » 

OIRAftPIG/DITJAOAT) 

RaRATE(JA-1)*DIP 
'10  KG«1 

lEfVLG.lT.VLECl))  GO  TO  93  ..  ... 

91  IECvUCKG)»Vt-G)92,93,94  _  _ 

- S7-Kl?iKe»l  . .  . 

CO  To  91 

93  OELGaPSUKG) 

-  GO  TO  90 

94  MIGA{VLG-VL£(XG-n)/tVLECKG).VlECKG-n) 

OELUA  CPSL^KG)'»0SL(EG«1))ahIG*P5L{KG«1) 

— 99'COimWE  . . 

lE(BP-CHG)  11,80,80 
8«  ONOT  30 
GO  TO  12 

11  IE  tIPT.NE.33  GO  TO  140 
lE(RT-OAN)  145,145,150 

143  R»i'R/l,S)»  XT,  ♦  C ,  5  •  R  T  /O  A  N)  7  — 

RT»R.T*RaOELTA 

130  :?(EPO-.9)  155,135,140 
159  3A»EA8ARI«n.-(EPUAAl.7n 
140  ONOToRaRMOaBARE 

9PaSP*(ONOTaDELTA) 

- I?  EHCR-p  ,tt;  cPfrsy'EfCPt5  'wm 

IE ( JEG.GT , 1 )  GO  TO  12 
IE CCPRS.-T.SABPR)  GO  TO  54 
OEETAa, 000005 
JEG*8 

12  V£L1»VEL2 

—  "trpr  rvrL31 - -  ■  ' 

RTM-CVEi,2*VEL2)/(GAMA.32.17.EIHP.PlT*TF/XTF) 

ROOTa1./(1.*CC(GA«A-1.)/2,)aBTMJ) 

TUTaROOT**TAD 

PREXaCPRSaTUT 
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TABLE  5.  CONTINUED 


<0*1 

so  If  (VEE  (K(1)-vEL2)51 , 52,53 
31  <0*<0»1 
GO  TO  30 

5?  PFACePST  CK03  ' 

GO  TO  59 

53  HAG*fVEL2-VEE(KO-13 3 / (VEE (K03 -VEE (KO-1 ) 3 
PFAC*  (  (PSY(K03  -PST  CKO-l}3*HAG)+PSY(KO-i3 
39  EFM*(SH0T*(ChG/PFaCJ)/32,17 
PRO**PRE* 

If CvEl2.lT.0.3  Prox*CPrS 

ACEL*  (PR0X-CPLG3*AREA*32.17/SH0T 
v£L2»vEH*(ACEL*0ELTA} 

ACOt*(ACEt*0ELTA)/2. 

3INC*(VEU*DELTA}*CACOL*DELTA3 

B0IS«B0IS*(BINC*12.3 

nr  AvoL.cvoLifBDis.AftEA)  —  .  * 

If { I LF -33  22,22,20 

22  If (RfST-CPLG)  19,19,23 

19  lLf-3 

20  »IGT»CPEG*32.17*»LEA/HGS* 

VLGI *VLG 

VLG»VLG»(AL5T. DELTA',  . . . ~  . — 

BLIC*(VLG.OELTAU(ALGT*OElTA*OELTA/2,3 

BU3«BLI3*BLIC*12. 

23  TLGS*  TLGS*OTDT  .  DELTA 
18  OLOT *VEL2 

glex*-gmlg/glg 

~  PLG8*l.*(GLG*VLG*VLG*HHSU7T84.3iTCG9iRTF*GHLG))  . . . 

PLGR*PLG0»*GLEX 

cpeg*cplg*plgr 

VLGS«HGIV*(BLI3.AuFA)-(8D1S.AREa) 

VLGS*XLGS-CVL*HG</3. 

If (OFLG.EO.0,3  PHG<*0. 

IF (OFLG,Ta.0,T  GO  TO  28  . . . . 

PMGP*HGH/(0FLG*32,  17} 

28  CONTINUE 

PLGa*12,»CPLG«aREa*VEl2/VLGS 

EFGL«(PHGM*{HGSH/32,17)3»CBLTA+1.3 

PLGB*12t*CPLG«AcEi«yLG/VLGS 

- P  L  G  C*T  21  iTrOLV VLTTiTLG  T»Et  G  /  V  L  G  5  . .  * 

PLGO*12.*MGM«RTF*DTOT/(GLG*YLGS«HMAL3 

DPLG*PLGA-PLGB-PLGC*PLGD 

™  ■  DTDT*  TLGS  •  OPtG  ‘•GAG'/CPl.G-  .  . . .  . . 

IF-(  ILF-3  J  54,54,56 

_ 5_4  XuE*«UE*l  _ _ _ _ _ _ _ 

GO  TO  37 
85  <U£»0 

“36  Continue 

<UV*KUV*1 

IFCKUV.LT. 2)  GO  TO  57 

- KUV.ff - - ~ 

WRITE (6,99  3  TIME , CPR3 , BO  I S , FPU, OPOT ,OLDT , CPLG, VLG, BLIS, CPEG , TLGS 
99  FORMAT (f 12.7,F12.1,f 12.3,f 12.4.7F1 1 .2) 

“57  CONTINUE  .  .  ' 

IF  (HLGP.LT.CPLG3  HLGP*CPLG 
CPLG«CPLG*(0PLG«0ELTA3 

- 0UDT“*9EL5 - 

UBW»CCHG-8P)/RH0 

COVL*CVL*BP 

GAMV«GAM0 

IFtACEL.LT. 0.3  GAMV.0. 


35 


TABLE  5.  CONCLUDED 


0PDT«( (0n0T«FI«P*12.)-( (G am v »E Fm» aCEL*VEL2» 1 2 . ) ♦ ( ARE A.CPRS* VEL 1  ♦  1 2 
1.)))/(Av0L-CU8**C0VL)) 

dope  * J 

TINEbTIME+OELTA 

Fey-BP/CRG .  " 

IF (CPRS.GT.PTOP)  PTOPbCPRS 
GN«(GINbFVOL*BP)  /  (AVOL-(U8w*COVL) ) 

TF«CPfi5«RN0L/(GN*RTF*lZ.5  . .  . . . . . 

CPRS»CPRS« (OPOT*OELTAJ 
IF(CPLG.GE.B).)  GO  TO  39 

HGBL.BLIS . .  . . ■ - 

WRITE(fl,4l)  9LIS 

41  F0RMAT(//47M  CASE  TERMINATED  AS  MAXIMUM  VELOCITY  REACHED  AT,F6,t, 

U7W  inches  of  travel///) 

39  CONTINUE 
79  CONTINUE 

'  CNITifBLT  S-hGBD/hGBL - - - - 

CPICb(vLG»VLGI)bCnIT*VLG 
WRITE (6,111)  CPIC 

MI  FORMAT(///10X,33H  LIGHT  GAS  GUN  MUZZLE  VELOCITY  IS, F8. 1 , 6HFT/SEC) 
WRITE (6,  1 12  3  MLGP 

117  FORMAT(lflX,28MLIGHT  GAS  PEAK  PRESSURE  I8,F8,1,4H  PSD 

- 8FTTETB;if3r  pcpp  - - - - 

113  FORMAT (IPX, 32HCGMPUSTI0N  SIOE  PEAK  PRESSURE  IS,F8,1,4H  PSI) 

GO  TO  « 

END  '  "  - .  . . 
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INITIAL  DISTRIBUTION 


I  IQ  USAF  (SAMID)  1 
HQ  USAF/SAMI  I 
AT SC  (DLTW)  1 
AT  SC  (SDl'A!  j  1 
ASD  iLNYS)  1 
FTD  (I'DXA)  1 
ATML  (LNP)  1 
TAC  (DRFM)  1 
TAC  (DM AM)  1 
WRAMA  (MMLBL J  1 
AliDC  (ARO.INC)  1 
AUL  (AUL- LLS-70- 239)  1 
CI1  ROD,  ATTN:  CRD  AM  1 
REDSTONE  SCI  INFO  CTR  2 
USA  KILN'S  COMO  A'lTN :  AMSWEREW  1 
ARMY  MATRL  SYS  ANALYS  AGCY,  AMXRD-AA  I 
USA  ABERDEEN  ROD  CTR,  AMXRD-BTL  2 
FRANKFORD  ARS,  ATTN:  LIB  1 
PICATINNY  ARS  AITN:  SMUPA-RT-S  1 
NAV  AIR  SYS  CMD,  CODE  AIR  3S0B  2 
USN  IVPNS  LAB  CODE  TR  1 
USN  ORD  LAB  (CODE  730)  2 
NAV  ORD  STN,  TECH  LIB  1 
NAV  SYS  CEN,  NEWPRT  LAB  1 
USNKC  (CODE  753)  2 
USNWC  (CODE  4585)  1 
USN  WPNS  EVAL  FAC  CODE  WE  1 
57  FTR  KPN  KG  (FWOA)  1 
OFF  CH  NAV  OPNS  1 
USN  RSCH  LAB  CODE  2027  1 
LOS  ALAMOS  SCI  LAB  1 
8ATTELLE  MEM  INST  1 
INST  FOR  DEF  ANALYSES  1 
SAND I A  CORP  2 
RAND  CORP  I 
USAFTFWC  (CRDC)  1 
1  LARRY  DIAMOND  LAB  1 
DDC  2 
TAWC-AY  1 
DL  1 
DLB  1 
DLOSL  2 
DLY  2 
DLR  1 
WE  1 
ASD  (YFM)  10 
DLDC  ID 
DLDL  30 
TRADOC/ADTC/DO  1 
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